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Abstract. We investigate the orbital selective Mott transition in two-band 
Hubbard models by means of the Gutzwiller variational theory. In particular, 
wc study the influence of a finite local hybridisation between electrons in different 
orbitals on the metal-insulator transition. 
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1. Introduction 

Metal-insulator transitions in Hubbard models with different densities of states have 
attracted a lot of interest in recent years [1-10]. A dispute arose over the question 
whether or not the transition occurs at different interaction strengths for the wide and 
the narrow band. A transition with different critical interaction parameters is usually 
denoted as an 'orbital selective Mott transition' (OSMT). Apparently, a consensus 
has been reached that such an OSMT can occur in Hubbard models with different 
bandwidths, subject to the bandwidth ratio a of the narrow and the wide band and 
the value of the local exchange interaction J. 

In most of the calculations in [1-10] the dynamical mean-field theory has been 
employed. We will use multiband Gutzwiller wave functions in order to study 
the OSMT. Such wave functions were originally introduced by Gutzwiller |llj in 
order to study ferromagnetism in the one-band Hubbard model. The evaluation of 
expectation values for the Gutzwiller wave function poses a difficult many-particle 
problem. Therefore, Gutzwiller, in his original work, used an approximation based 
on quasi-classical counting arguments [12l[13]. This 'Gutzwiller approximation' later 
turned out to be equivalent to an exact evaluation of expectation values in the limit 
of infinite spatial dimension or infinite coordination number. Generalised Gutzwiller 
wave function for multi-band Hubbard models have first been introduced and evaluated 
in the limit of infinite spacial dimensions in reference [15| . The formalism was further 
generalised, e.g., for superconducting systems, in references [TBI [TTj . 

The OSMT in a two-band Hubbard model has first been investigated by means 
of the Gutzwiller theory in reference [9]. In that work the authors found an OSMT 
both for vanishing (J = 0) as well as for finite (J ^ 0) local exchange interaction. 
For J = the critical band width ratio was found to be ac ~ 0.2. The Gutzwiller 
results in [3] were in good agreement with data from DMFT and a slave-spin approach 
proposed in reference [7]. 
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In this work wc will analyse the OSMT in a two-band model in more detail. 
In particular, we permit a finite expectation value Ag = (c|^]^Ci.2) for the local 
hybridisation which can change the nature of the OSMT. Such a hybridisation could 
be finite spontaneously, solely due to the Coulomb interaction, or due to a finite 
hybridisation term in the Hamiltonian. We will investigate both possibilities. 

Our paper is organised as follows: The two-band Hubbard models are introduced 
in section [21 In section [3] wc define generalised Gutzwiller wave functions and give the 
results for the variational ground-state energy for these wave functions in the limit of 
infinite spatial dimensions. The orbital selective Mott transition in a two-band model 
without a finite local hybridisation is discussed numerically, and as far as possible 
analytically, in section [4l In section [5] we investigate analytically the spontaneous 
hybridisation in a spinless two-band model. Finally, the hybridisation effects in the 
full two-band are studied in section [6l and a summary closes our presentation in 
section [71 

2. Model systems 

In this work we investigate the two-band Hubbard model 

^ = E 4ms-,^<. + E = ^0 + ■ (1) 

i.j:b:(7 i 

Here, the one particle Hamiltonian Hq describes the hopping of electrons with spin a 
on a lattice with L sites. The index 6=1,2 labels the two degenerate orbitals at each 
lattice site. We assume that the hopping amplitudes 

^tj = C'bUj (2) 

depend on the orbital index b only via overall bandwidth factors at- This leads to an 
orbital-dependent renormalisation 

D,ie) = -Do (-) (3) 
ah \ctb J 

of the bare density of states 

i?o(e) = ^E'^(^-^fc)' (4) 

k 

where Sk is the Fourier-transform of tij. Throughout this work, only symmetric 
densities of states will be considered Do{—e) = Do{s). 

We will study the two-band model ([T]) with and without spin-degrees of freedom. 
For the full two-band model we assume that the orbitals have an eg-symmetry. The 
atomic Hamiltonian then reads 

i?l?^ = c/E 

«6,T"-fca + f^' E "1,'^"2,<t' - -'^ E "1.'^"2,<t (5) 
b a.a' cr 

a 

where in cubic symmetry the two parameters U' and Jc are determined by t/' = U—2J 
and Jc = J- Without spin, the atomic Hamiltonian -ffi:at simply reads 

hH^ = Uh,n2 , (6) 
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where the effective Hubbard interaction in this model can be derived from the 
interorbital Coulomb {U') and exchange (J) interaction through U — U' — J. 
Apparently, the spinless two-band model is mathematically equivalent to a one-band 
model with a spin-dependent density of states. In the limit q;2 ^ it becomes a 
Falicov-Kimball model. 

Both atomic Hamiltonians @ and ([6]) can be readily diagonalised 

i7«'(')=^i?r|r)(r|. (7) 
r 

The eigenstates |r) of iJ^j-* are the empty state |0), the two singly occupied states \b) 

^ (2) 

and the doubly occupied state \d). The diagonalisation of H)^^ leads to similar Slater- 
determinants for all particle numbers riat ^ 2. In the two-particle sector, riat = 2, one 
finds the triplet ground-state with energy Er = U — 3J, in agreement with Hund's 
first rule, and three singlet states with energies Er = U — J (doubly degenerate) and 
Er = U + J; for more details, see reference [15]. 



3. Gutzwiller wave functions 



3.1. Definition 

In order to study the two-band Hubbard models introduced in section [51 we use 
Gutzwiller variational wave functions [llj which are defined as 

1^0)^11^0) ■ (8) 

i 

Here, |5'o) is a normalised one-particle wave function and the local correlation operator 
Pi has the form 

P = ^ Xr,r''rrirx' , (9) 
r.r' 

for each lattice site i, and 

mr,r' = |r)(r'| . (10) 

The real coefficients Ar.r' and the one-particle wave function |4'o) are variational 
parameters. For systems without superconductivity it is safe to assume that the 
parameters Ar,r' are finite only for atomic states |r), |r') with the same particle 
number. For ground states without spin order one can further assume that only states 
with the same Sz quantum number lead to finite non-diagonal variational parameters. 
Due to these symmetries the correlation operator ^ contains up to 5 variational 
parameters for H^l^ and up to 26 for //^t''- 

Throughout this work we will investigate the half-filled case of our model systems 
and allow for a finite local hybridisation 

^0 = (cl,l,rTC,;,2,^>*0- (11) 

With respect to the operators and c, the local density matrix is therefore non- 
diagonal. For analytical and numerical calculations, it is more convenient to work 
with creation and annihilation operators 

W) 1 /'-At) , i(t) A 
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which have a diagonal local density matrix, 

nH"^ = {hlb.Ab'J^o = SbM ± Ao) . (14) 
With these operators the one-particle Hamiltonian Hq reads 

^= E ^ihl^Ab'.. (15) 

i,j;b,b';(T 

where 

i--'' ^^-fiSb.b' +Aail-5t^i,')) ■ (16) 

Both atomic Hamiltonians ^ and © keep their form under a transformation from 
c to h. By building a basis of Slater determinants \H) with the operators /i- f,^ the 
eigenstates of the atomic Hamiltonian can be written as 

|r) = 5]rr,H|i?). (17) 

H 

3.2. Evaluation in infinite spatial dimensions 

The evaluation of expectation values for Gutzwiller wave functions poses a difficult 
many-particle problem. In this work we employ an evaluation scheme that becomes 
exact in the limit of infinite spatial dimensions. Within this approach the expectation 
value of the local Hamiltonian reads 

(^at)*G = E -^rAri,rAr,r2("T-ri,r2)'i'o • (18) 

Here, the expectation value ("^ri.r2)*o given as 

(TOri,r2)*o = E^ri,_f/Tr2,ifm^ (19) 

H 

where 

<= n "i'^ n (i-"i'^)- (20) 

b{occ.) 6(unocc.) 

For the expectation value of a hopping term in the one-particle Hamiltonian one finds 

('iI,b,a'lj,b%<T)*G '='^%b%'b'iHAa'^],b',a)'S'o , (21) 
b,b' 

where the elements of the renormalisation matrix q are given as 

\^jr4)(r 



Ibb- E Ar,x2Ar3,r.(r.|fet,jr3) ^^ ' (22) 



1 — 



The remaining expectation value in ((22)) can be calculated in the same way as (|T9)) . 
Note the symmetries qi i — ^2,2 and gi,2 — 92,1- The renormalisation factors for the 
c-operators are diagonal, 

(^i,f),(TS,fc,<^)*G — ^b{'^\,b,a^j,b,a)^o (23) 

and given by 

9(1-) = 91,1 ± 91,2 ■ (24) 
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Furthermore, the evaluation in infinite dimensions shows that the variational 
parameters Ar,r' a-nd the one-particle wave function jV^o) have to obey the constraints 

l = (p2)^,„= ^ Ari.rAr,r2("iri,r2)*o (25) 
r,ri,r2 

and 

nt'^S,.,, = {P%/KJ^„ (26) 

= X! ^ri,rAr,r2(™ri,r2^6.cr^b',cr)*o ■ 
T.V1.V2 

4. The orbital selective Mott transition in a two-band Hubbard model 

In this section we investigate the metal-insulator transition in the two-band Hubbard 
model without local hybridisation. We use a semi-elliptic density of states 

D^{e) = -Vl-e^ (27) 

TT 

which leads to the bare one-particle energy 

f" 2 
£0= / deDo{e)e^-— . (28) 

Our energy unit is given by I? = 1, half of the bare bandwidth. When we set ai = 1 
and introduce the bandwidth ratio a = ai/a2 < 1, the expectation value for the 
one-particle Hamiltonian in ([T]) is given as 

(iJo>*G = (9?+'?2a)£o. (29) 

Without hybridisation, the variational ground-state energy has to be minimised only 
with respect to the variational parameters Ar.r'- In figure[l](left) we show the resulting 
renormalisation factors qb as a function of U for J = and two different bandwidth 
ratios a. As already observed in reference [5], it depends on the value of a whether 
or not there is an orbital selective Mott transition. For J = 0, the critical ratio is 
ac = 0.2, i.e., the renormalisation factors gi, q2 vanish at two different critical values 
Uc2 < Uci if a < etc- By switching on J, the critical ratio ac becomes larger and the 
Mott transitions take place at smaller values of U ; see figure [TJright). 

For J = 0, we can gain more insight into the nature of the different Mott- 
transitions in our model by some analytical calculations. First, we consider the 
case a > ac- If we approach the Mott transition from below, we can neglect the 
variational parameters niQ = rui for empty and fourfold occupied sites. Due to the 
high symmetry of the model for J = the ground-state energy is then a function of 
only three variational parameters, d, </>, and, 6*, 

E = 2eod{l-2d)f{(l),e) + {l + d)U (30) 

where 

/(0,6') =4ai (sin((/))sin(e') + ^/2eos((/))cos(6'))^ (31) 

2 



4a2 (cos(0)sin(6') + V2sin((/))cos(6l)) . (32) 



Here, tan((/))^ gives the ratio of the probabilities to find a singly occupied site with 
an electron in the wide and in the narrow orbital. The ratio of the probabilities 
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for doubly occupied sites with two electrons in the same and in different orbitals is 
parametrized by tan (0)^. The variational parameter d gives the total probability for 
single occupation. At the Mott transition, where c? ^ 0, the two angles 0, 9 can be 
calculated analytically 

6o = e{d 0) 
(/)o = <t>{d ^ 0) 

Both values, tan ((/)o)^, and tan (^o)"^ are shown as a function of a in figure [2jleft). 



1 /-17 + 2a-17a2\ 

iarctanfii±^l^^^^^) f34) 
2^'^^*^'^((l_c.)(l + cos20o)^' ^ ' 
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Figure 2. left: tan (0o)^ (straight) and tan (9o)'^ (dashed) at the Mott transition 
as a function of bandwidth ratio a; right: Ratio of renormalisation factors (72/91 
at the Mott transition as a function of a. 



As expected, the weight of local states with no electron in the narrow band 
vanishes for a ^ a^. The renormalisation factors qi, both vanish proportional to a 
V^CTfT^^TJ, when U approaches Uc from below. The ratio (?2/'?i is 

see figure [2] (right). 



square-root, 
finite for U 



Uc and goes to zero proportional to ^Ja — ac 



Finally, the critical interaction strength Uc2 = Ud is given as 
Uci = 2|eo|/(0o,6'o) (a > a^) . 



(35) 
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Next, we consider the case a < a^. For interaction parameters 11^2 < U < Ud, 
the electrons in the narrow band are locahsed and the wide band can be treated as an 
effective one-band model. This leads us to the critical interaction parameter 

[/ci = 2|£o|/(0, 0) = 16|eo| (a < «c) (36) 

for the Brinkmann-Rice transition of the wide band. Starting from the Brinkmann- 
Rice solution for U < Ud, we can expand the variational energy to leading (i.e. 
second) order with respect to the three parameters {vi} = {(j), 9, mg}, 



E = Eo 



3 

J2 v^Ei^jVj . 



(37) 



The localisation of the narrow band becomes unstable when the matrix E has negative 
eigenvalues for physical parameters Vi > 0. This evaluation yields the following 
expression for the narrow-band critical interaction strength 

C/c2 = 16|e|-^ (a<ac). (38) 
1 — 4a 

The resulting phase diagram for all < a < 1 is shown in figure [3] 




0.4 0.6 
a 



Figure 3. Critical interaction parameters Ud (straight) and Uc2 (dashed) as a 
function of a (see eqs. I|35|l. I|36|l. l|38})- 



5. The spinless two-band model 



As the simplest example for a model with different densities of states we investigate the 
spinless two-band model. In the half filled case and without spontaneous hybridisation 
(Aq = 0) the constraints ipS]) and (|26p can be solved analytically for this model. The 
variational energy is then solely a function of A^, 



E^ 



= 4A^ ( 1 - 



A5 



£o + -tA, 



(39) 



The energy p9p can be minimised analytically. As a result one finds the well known 
Brinkmann-Rice solution 

.BR, = l-(^)', (40) 
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4 V C/c, 

for the renormalisation factor q = Sb^b'Qb.b' and the expectation vahie of the double 
occupancy d = A^/4. The Brinkmann-Ricc metal insulator transition occurs at the 
critical value U = Uc = 16|eo|- 

For the renormalisation factors we set ai + a2 = 2, i.e. the difference of the 
bandwidths is parametrized by Aa = a\ — a-i. Starting from the analytic solution 
for vanishing hybridisation we can calculate the variational ground state energy to 
leading order in Aq, 

Sao = Sbr + C(C/, Aa)A2 . (42) 

A spontaneous hybridisation will appear if the coefficient C in (|42|) is negative. The 
analytical evaluation leads to the Stoner-type instability criterion 

/(Aq) [//£/a2 + _ 

UM^)^ 2{i + u/u.f ^^'^ 

where 

/(Aa) ^ . , ^ (44) 

2arcsmh (Aa/V4 - Aa^j 



0.5 




5 



Figure 4. left: /(A«); rigftt: g{U/Uc) 

In figure [5] the function /(Aa) and the right hand side of equation (|43|) are shown as 
a function of Aa and respectively. As can be seen from this figure the function 
/(Aa) and therefore the left hand side of (|43)) approach zero for Aa 2. On the 
other hand, the right hand side of (j43|) is positive for all U > 0. This means that for 
arbitrary values of U there exist a critical bandwidth difference Aac with Aq > for 
a > ac. Figured] (left) shows the phase-diagram for ground states with and without 
finite hybridisation for different values of the density of states -Do(O) at the Fermi- 
level. Whether or not there is a transition in the large U limit for all values of Aa 
depends on the value of -Do(O). This is illustrated in figure [5] (right) where the critical 
difference Aac for the transition is shown as a function of Dq{0) in the limit U ^ oo. 
Note that a spontaneous hybridisation has already been observed in a Falicov-Kimball 
model within a mean-field approximation [18| . This is in agreement with our results 
in the limit Aa 2. 

In summary, our analytical results on the spinless two-band Hubbard model 
show that a difference in the bandwidth increases the tendency of the system to 
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Figure 5. ie/t:phase diagram of the spinless two-band Hubbard model for 
different densities of states at the Fermi level D(0) = 0.25,0.2,0.15, 0.13,0.125, 0.1 
(from the bottom to the top of the figure) right: critical difference Aoc in the 
limit f/ — > oo as a function of Dq{0). 



exhibit spontaneous hybridisation between the narrow band and the wide band. 
Mathematically, the reason for this is quite simple. Both, the expectation value of the 
one-particle energy Tfg and the Coulomb interaction H\oc are changing quadratically 
in Aq. However, in the limit Aa the energy gain from H\oc always beats the rise 
in energy due to Hq. At first glance, one might think that the same behaviour should 
be observed in the OSMT phase of the two-band model with the only difference that 
it is not the bare but the effective width of the narrow band that vanishes. As we will 
discuss in the next section, however, this hypothesis turns out to be incorrect. 

6. Hybridisation in the two-band model 

In this section we present numerical results for the two-band model with a finite local 
hybridisation (jlip . The hybridisation can develop either spontaneously, like in the 
spinless model (section [S]) , or it can be caused by a finite hybridisation term in the 
Hamiltonian. We will discuss both effects separately. 

6.1. Spontaneous hybridisation 

As shown in section [SJ a vanishing width of the narrow band can be the driving force 
for a spontaneous local hybridisation of the wide and the narrow band. In our two- 
band model, however, the vanishing of the effective bandwidth for (72 ^ does not 
have the same effect. This can be seen in figure IHl where we show the results for the 
renormalisation factors 51, (72 and the hybridisation Aq. Unlike in the spinless model, 
there is not necessarily a finite hybridisation if the effective narrow bandwidth goes to 
zero for U —>■ Uc2- The reason for this differing behaviour is an additional contribution 
to the one-particle energy of the full two-band model. To leading order in Aq there is 
a third term from the expansion of the narrow-band renormalisation factor 

g2«'Z2(Ao = 0) + cAg. (45) 

The coefficient c is negative and, multiplied with the negative bare one-particle energy 
of the narrow band it leads to an increase of the total energy. This contribution to 
the energy overcompensates the negative term from the Coulomb interaction. 
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Figure 6. Rcnormalisation factors gi, q2 and hybridisation 2Ao for J = and 
a = 0.15. 



A finite hybridisation Aq sets in at larger values of U when the system is already in 
the OSMT phase, see figure[6l Numerically, it seems as if Aq approaches its maximum 
value A^'^'' ==1/2 only in the hmit U ^ oo. 

In all systems with finite values of J that we investigated, we did not find a 
solution with spontaneous hybridisation. It is more likely, though, that for values of J 
smaller than some critical parameter there is a solution with a finite hybridisation. 
However, it is difficult to determine this small parameter numerically. 

6.2. Finite hybridisation in the Hamiltonian 

The assumption that there is no hybridisation between the two degenerate bands in 
the Hamiltonian of our model is quite artificial. In this section we will therefore 
investigate how the OSMT is affected if we add a hybridisation term of the form 



to our Hamiltonian For J = we find that the OSMT phase is destroyed for any 
finite value of fj. This is illustrated in figure [7] (left) where we show the expectation 
value Aq as a function of U for several values of fj. 

For finite J, the behaviour of our model is more ambiguous. As we have seen 
before, a finite J stabilizes the OSMT phase whereas a finite f] tends to destroy it. 
Therefore, it depends on the ratio of both quantities whether or not an OSMT is 
found. Figure [7] (right) shows the renormalisation factors qi, for different values of J 
and 77. For J = 0.025C/ and f) = 0.05D the OSMT is completely suppressed. This 
is still the case for the smaller value fj ~ Q.025D, although the narrow band factor 
q2 is already quite small in the region of U parameters where it would be zero for 
77 = 0. Finally, for larger values J = 0.05U an OSMT phase is restored for interaction 
parameters U > Uc2 where Uc2 is larger then the corresponding value for fj = Q. 

In summary, our numerical calculations show that appearance and disappearance 
of an OSMT results from a subtle interplay of the local exchange interaction J and 
the local hybridisation i). 




(46) 
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Figure 7. ie/t:expoctation value 2Ao as a function of U for several values of i); 
rig/it:Renormalisation factors qi, q2 for a = 0.15 and jj = 0.025D, J = 0.05U 
(straiglit), ij = 0.025Z), J = 0.025(7 (daslied), i) = 0.05D, J = 0.025{/ (dotted). 



7. Summary 



In this work wc have investigated the orbital selective Mott transition (OSMT) in two- 
band Hubbard models with different densities of states by means of the Gutzwiller 
variational theory. We were particularly interested in the question how the OSMT 
is modified when we allow for a finite local hybridisation between the wide band 
and the narrow band. In the two-band model without spin-degrees of freedom there 
always is a spontaneous hybridisation if the narrow bandwidth goes to zero. However, 
we did not find such a behaviour in the full two-band model. There, spontaneous 
hybridisation was only seen for vanishing local exchange interaction, J = 0, and for 
Coulomb parameters U larger then the critical parameter at which the electrons in the 
narrow band localise. By adding a local hybridisation term ^ fj to the Hamiltonian, 
the phase diagram becomes more involved. Whether or not an OSMT takes place 
depends on the relative strength of J and fj. The exchange interaction J tends to 
stabihse the OSMT phase, whereas the hybridisation 77 tends to destroy it. 
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